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Mineral inclusions in biomass are attracting increased scrutiny due to their potential impact on processing 
methods designed to provide renewable feedstocks for the production of chemicals and fuels. These 
inclusions are often sculpted by the plant into shapes required to support functional roles that include the 
storage of specific elements, strengthening of the plant structure, and providing a defense against pathogens 
and herbivores. In situ characterization of these inclusions faces substantial challenges since they are 
embedded in an opaque, complex polymeric matrix. Here we describe the use of Bragg coherent diffraction 
imaging (BCDI) to study mineral inclusions within intact maize stalks. Three-dimensional BCDI data sets 
were collected and used to reconstruct images of mineral inclusions at 50-100 nm resolution. Asymmetries 
in the intensity distributions around the Bragg peaks provided detailed information about the deformation 
fields within these crystal particles revealing lattice defects that result in distinct internal crystal domains. 

Characterization of a heterogeneous sample such as corn stover involves dealing with material that is 
organized on multiple length scales and contains discrete domains with different compositions, orienta- 
tions and order. Investigations into the structural components of these types of sample are of considerable 
interest since the presence of mineral inclusions in plant tissues is a factor in the development of industrial 
processes used to produce chemicals and fuel 1 . In addition, the inclusion of crystals in plant tissues has been 
identified with several distinct biological functions that include providing sinks for chemical accumulation, 
adding structural support to plant tissues and creating a defense against pathogens and herbivores 2-4 . 

A novel approach for non-destructive in situ imaging of the sizes and morphologies of micro crystals embedded 
in opaque media is provided by Bragg coherent diffraction imaging (BCDI). This approach has previously been 
used to characterize manufactured samples of simple crystals 5 " 8 . Here we demonstrate the applicability of BCDI to 
more complex crystals embedded in biological samples and, more generally, the use of the technique for analyzing 
sample types in which crystalline particles are concealed within an amorphous media. 

Results 

In order to survey the nature of x-ray scattering from the vascular bundles embedded in maize stalks, we collected 
wide angle x-ray scattering (WAXS) data using an 180 mm sample-to-detector distance at the BioCAT beam line 
at the Advanced Photon Source (APS). These patterns exhibited the characteristic intensity distribution expected 
for cellulose ip 9 as well as occasional, very sharp reflections of uncertain origin (Fig. 1). These sharp reflections 
exhibited no apparent preferred orientation relative to the fiber axis of the stover. Thirty WAXS patterns from 
maize stalks were examined and the scattering angles of nearly 100 of these Bragg peaks were measured from these 
data sets. A comparison of the measured scattering angles with the d-spacings of reflections predicted from the 
crystal structures of naturally occurring minerals commonly observed in plant tissues concluded that almost all of 
these reflections could be indexed using the crystal lattices of calcium oxalate monohydrate (whewellite) and 
dihydrate (weddelite). The strongest reflections expected 1011 from whewellite (5.93, 3.64 and 2.98 A) and wed- 
delite (6.19, 4.42, 3.68 and 2.75 A) were all observed in this survey, indicating that both crystal forms exist in these 
samples. The observation of Bragg peaks on most of these images shows that these microcrystals occur at a density 
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Figure 1 | WAXS Pattern from vascular bundles of maize. This WAXS 
pattern shows diffraction characteristic of cellulose fibers and containing 
several Bragg reflections (magnified in the insets) that are due to 
diffraction from crystallites embedded in the material. At the point-to- 
point resolution of this pattern, taken at a sample to detector distance of 
180 mm, little detail within the Bragg reflections can be discerned. 

that is typically high enough to find crystals in orientations that place 
a few strong reflections on the detector surface. For the x-ray beam 
size of 140 X 40 |im and assuming a sample thickness of ~ 1 mm the 
sample volume illuminated by x-rays in these experiments is ~6 X 
10 6 |im 3 . The calculations of three-dimensional diffraction patterns 
from the crystal structures 10 of whewellite and weddelite showed that 
they produce very few strong reflections and, therefore, the probabil- 
ity of detecting diffraction from a specific randomly oriented crystal 
is small. If the mosaic spread of a crystal is assumed to be —0.05° the 
crystal number density probably must exceed one crystal per 10 4 |im 3 
sample volume to produce detectable Bragg peaks on the majority of 
WAXS images. 

X-ray fluorescence microscopy (XFM) revealed occasional cal- 
cium-rich inclusions with linear dimensions of ~1 micron (Fig. 2). 
Since the samples used for the XFM analysis were from the same 
source as the stover used for the WAXS analysis these observations 
are consistent with the presence of calcium oxalate crystals. Although 
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Figure 2 | XFM images of maize. XFM images of a fiber cell teased from a 
sample of corn stover. Cell walls around the fiber cells are typically 10 
microns in diameter and hollow. The distribution of most elements 
corresponds to nearly uniform distribution within the cell walls. The 
exceptions are calcium which exhibits several micron- sized calcium-rich 
inclusions, and silicon which appears to concentrate within the hollow 
interior of the fiber cells. 



the resolution of the XFM images is relatively low, the inclusions 
appear to have well-defined shapes and edges. Silicon-rich inclusions 
were also observed, particularly in the hollow centers of fiber cells. 
However, Bragg peaks observed by WAXS could not be indexed on 
the basis of any known naturally occurring silicon-rich mineral. We 
concluded that the silicon-rich inclusions were either amorphous or 
constituted crystals sufficiently small that the scattering peaks were 
not detectable in the WAXS patterns. 

BCDI diffraction data were collected on a CCD detector mounted 
on a 6-circle diffractometer at beam line 34-ID-C at the Advanced 
Photon Source at Argonne National Laboratory using a beam size of 
approximately 2 (im 6,12 . Early observations at room temperature 
indicated that radiation damage to the surrounding organic material 
was leading to movement of the crystals in the beam. Consequently, 
all data were collected with the samples held at liquid nitrogen tem- 
perature. Vascular bundles teased from dried maize stalks were 
adhered to the cold finger of a cryostage with their fiber axes per- 
pendicular to the incident x-ray beam. The individual vascular bun- 
dles, typically 1-3 mm in diameter, were broken from the bulk of the 
stalk and laid parallel to each other on the copper surface of the cold 
finger (see Supplementary Information, SI). The detector was placed 
at a scattering angle corresponding to a spacing of 3.65 A where 
scattering from both forms of calcium oxalate can be observed. 
The samples were translated across the x-ray beam with the detector 
at 0.75 m from the sample in order to provide a relatively large field 
of view for detection of Bragg reflections. Once a Bragg reflection was 
observed, the detector was typically translated to either 1.0 m or 
1.8 m from the sample in order to collect the data on the finest 
possible grid. By rotating the sample over small increments a 




Figure 3 | BCDI data, (a) Sections through a three-dimensional BCDI 
data set collected as a series of diffraction patterns in which each successive 
pattern is separated from the previous by a small, constant angular rotation 
of the detector, 56. This illustration shows every fifth image taken as the 
detector rotated through the reflection. The stack of almost parallel images 
is then assembled into a three-dimensional data set that maps the detailed 
diffracted intensity around the position of a single Bragg reflection, 
(b) Central sections through five different Bragg peaks corresponding to 
the five reconstructions shown in Figure 4. The intensity measurement in 
all of the diffraction images are displayed on a logarithmic scale. 
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three-dimensional data set was assembled as a set of two-dimen- 
sional slices that map the intensity through the reflection (Fig. 3a). 

About a dozen data sets suitable for three-dimensional reconstruc- 
tion were collected from vascular bundles harvested from dried 
maize stalks. Figure 3b shows a gallery of central sections through 
the diffuse speckles around Bragg peaks from five different samples. 
Measuring the distance from the center of a reflection to positions 
where the intensity disappears into the background suggests that the 
resolution achieved by the most strongly diffracting samples was 
—50 nm in the dimensions that correspond to the detector surface, 
with more typical data sets achieving a resolution closer to — 100 nm. 
In the dimension corresponding to the sample rotation the max- 
imum resolution is constrained by the extent of the rotation and 
proved lower, typically in the range of —200-400 nm. 

Image reconstructions of calcium oxalate crystals utilized these 
BCDI data sets and a local implementation of algorithms previously 
demonstrated to accurately determine the sizes and shapes of crystal 
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Figure 4 | Crystal reconstructions from BCDI data, (a) Top and side 
views of isosurface renderings of the absolute-valued electron density for 
five crystals reconstructed from BCDI data. The threshold value of the 
surface contour was set to make visible the largest extraneous noise feature 
in the map. Noise and spurious detail in the reconstructed images were 
suppressed by local smoothing of the density over volumes with 
dimensions approximately equal to the resolution of the data. The inset 
image depicting the atomic structure of the whewellite crystal was obtained 
from the Mineralogy Database (http://webmineral.com/data/ 
Whewellite.shtml) 28 with the c-axis running in the vertical direction, 
(b) Orthogonal slices showing phase maps through the centers of the 
corresponding crystals shown in (a). The top row shows the set of cross- 
sections as viewed down the c-axis and the bottom row shows the two 
corresponding longitudinal sections for each reconstruction. The color 
table rotates from blue (0 radians), through red (n radians) and reverts to 
blue (271 radians). The values of phase angles outside the volume of the final 
support have been eliminated for clarity. The density isosurfaces were 
rendered using the Chimera program 29 and the phase diagrams were 
rendered with the Image] software 30 . 



objects 13 " 21 . The implementation was checked by reconstructing 
known objects from simulated data. The general approach involves 
the iteration of real- space object constraints, in which the image of 
the object under reconstruction is retained only within the confines 
of a support volume, and a reciprocal-space step, in which phase 
angles from the back transform of this modified image are recom- 
bined with the experimentally determined diffraction amplitudes. 

Renderings of crystals reconstructed from the BCDI data (Fig. 4a) 
indicate that they range in length from 2-5 |im with cross sections 
0.5-1 urn across. The images are oriented identically in the frame of 
reference of the detector with the plane of the detector horizontal and 
perpendicular to the plane of the images. The two longest crystals 
have a symmetric appearance with flattened cross sections and poin- 
ted ends whereas the other three examples are —2 um in length and 
have blunt ends. Some aspects of each crystal show indications of the 
presence of flat faces. Disordered material or crystalline material 
oriented in such a way that it does not scatter into the Bragg peak 
being used for the reconstruction will not contribute to the recon- 
struction and result in a void. Consequently, these images corre- 
spond to ordered crystalline domains giving rise to scattering in 
the observed Bragg reflection. 

The crystals, reconstructed on the basis of a Bragg peak at —3.65 A 
spacing, could be either calcium oxalate monohydrate or dihydrate 
since the (0 4 0) peak of whewellite and the (0 0 2) peak of weddelite 
both scatter very near that angle. However, an argument based on the 
known crystal symmetries for these two crystal forms shows that all 
of these reconstructions correspond to the monohydrate form. First, 
the 3.65 A periodicity is fixed in a direction perpendicular to the long 
axis of these crystals. In weddelite this d- spacing corresponds to the 
(0 0 2) reflection so the periodicity is placed on the crystal c-axis. 
Since weddelite crystallizes in a tetragonal space group the a- and b- 
cell axes are equal and the (1 0 0) and (0 1 0) directions are equivalent. 
Consequently, the required orientation would place one of these two 
equivalent axes along the long axis of the reconstructed crystal and 
the other perpendicular to the long axis. In contrast to this contra- 
dictory situation, whewellite crystals are frequently observed as hav- 
ing long, thin shapes with pointed ends. These crystals are referred to 
as 'raphides' and the c-axis of whewellite has been shown to corre- 
spond with the long axes of raphides 11 . These observations indicate 
that the crystals that we have reconstructed are whewellite. 

The x-ray scattering patterns about the Bragg peaks analyzed in 
this study are characterized by strikingly asymmetric intensity dis- 
tributions that sometimes contain a distinct substructure of streaks 
and an ill-defined central maximum (Fig. 3). Although the overlap- 
ping peaks that occur in merohedrally twinned crystals may also lead 
to complicated intensity distributions 22 , asymmetry in the intensity 
distribution around the position of a Bragg reflection generally 
requires that the real- space object be considered as a complex rather 
than purely real electron density function 8 ' 7 . A complex- valued elec- 
tron density may be represented by a function that contains both a 
magnitude and phase angle component. The magnitude of this func- 
tion represents the physical electron density and, as in the more 
conventional case of a purely real electron density function, the 
absolute value of this term defines the shape and extent of the object. 
The phase component measures the local deformation of the crystal 
lattice and variations in phase across the object correspond to shifts 
of diffracting units within the object relative to the periodicity of the 
ideal crystal. Consequently, a phase map calculated over the object 
volume may be used to identify positional variations in the crystal. A 
phase shift of n radians is indicative of a part of the object shifted by 
half a unit cell from its ideal position. For example, in BCDI analysis 
involving manufactured specimens of very simple crystals this rep- 
resentation has been used to map strain fields across the crystal 8,23 . 
Phase angles are sensitive to positional shifts on length scales less 
than the unit cell of the crystal and, therefore, probe changes in the 
crystal periodicity on an angstrom length scale. 
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Examination of the distribution of the phase angle component in 
the complex electron density reveals sharp separation of the crystal 
into domains of relatively uniform phase (Fig. 4b). This arrangement 
of domains is characterized by a phase shift exceeding nil radians, 
indicative of over a quarter-unit cell shift in one part of the crystal 
relative to another. Most of these reconstructions show a central 
crystal core surrounded by a shell with shifted phase. The presence 
of an intermediate structural layer of sub -resolution dimensions 
(—50 nm) at the domain boundaries cannot be ruled out although 
there is no evidence for a component of this type. AFM studies of 
whewellite crystal growth identify the crystal face associated with the 
crystal a-axis as the slowest growing face and prone to a variety of 
crystal imperfections 24 . The surface of this face is not atomically 
smooth since it is composed of oxalate ions in two different orienta- 
tions and a mismatch of crystal blocks at this interface (as in a grain 
boundary) would result in a significant phase shift across the crystal. 
In addition, when viewed down the b-axis the whewellite crystal 
contains planes of calcium and oxalate parallel to the c-axis in which 
a significant slippage would result in relatively little disruption in the 
organization of intermolecular bonds (inset in Fig. 4a). 

The x-ray scattering at a particular Bragg reflection (in this case the 
(0 4 0) reflection) is associated with a specific orientation of the 
reflecting planes in the crystal relative to the incident x-ray beam. 
For this reason the lattice orientation in the crystal reconstructions 
from the BCDI data is fixed by the selection of Bragg reflection for 
analysis. For the largest crystal in Fig. 4a the c-axis is vertical although 
this orientation is tilted by about 20° from the long axis of the crystal. 
The phase distribution in this crystal (Fig. 4b) indicates that the 
dislocation plane is parallel to the c-axis of the crystal lattice rather 
than the long axis of the crystal. This suggests the dislocation is not 
uniquely linked to the crystal habit but is defined by the orientation 
of the underlying lattice. 

Discussion 

Plants control the size and shape of biomineral inclusions 25 . 
Although capable of taking on several distinct habits, whewellite 
most often takes on the elongated shape known as a raphide. 
Raphides in Araceae are elongated, frequently pointed, and accom- 
modate grooves along their long axes 3 . The combination of sharp 
ends and adsorbed proteases produces toxic responses when plants 
harboring raphides are ingested. The sharp ends of these crystals 
appear to be developed so as puncture cells and introduce irritant 
into organisms ingesting the material 2 ' 3 . The observation here of 
crystal particles organized into distinct crystal domains related by 
a large dislocation of the crystal unit cell is suggestive of a new level of 
biological control in the deposition of biominerals. The function of 
this dislocation is not clear but may contribute to the formation of the 
sharp ends of raphides. 

Cellulosic ethanol production technologies make use of the whole 
plant (rather than just the corn kernels) and are, therefore, poten- 
tially appealing to the biofuel industry as a more efficient means for 
generating ethanol than standard fermentation techniques. Calcium 
oxalate crystals are the prevalent mineral deposit in higher plants but 
have usually been identified in plant leaves and have not, to our 
knowledge, been previously reported in maize stalks. Characteri- 
zation of the chemical composition and sizes of calcium oxalate 
particles embedded within maize stalks could be of value in the 
development of industrial processes for their breakdown and sub- 
sequent purification to ethanol. Additional techniques, such as x-ray 
micro -tomography, are able to provide the number density of micro 
crystals in plant tissues 26 . 

These results demonstrate the level of detailed characterization that 
can be accomplished in the study of microcrystals embedded in an 
opaque, heterogeneous material; identification of the mineral com- 
position of the inclusions; their size and shape; the orientation of the 
crystal lattice within the crystal; the form and nature of dislocations. 



The microscale architecture of heterogeneous materials, natural or 
synthetic, is central to providing the extraordinary properties that 
make them increasing popular for a broad range of applications. 
BCDI represents a novel, non-invasive approach to characterization 
of this architecture at a level of detail difficult to achieve by other 
methods and essential as a basis for future materials design. 

Methods 

WAXS data were collected at beam lines 18-ID (BioCAT) and 23-ID-B (GMCA CAT) 
at the Advanced Photon Source (APS), Argonne National Laboratory, using a sample 
to detector distance of approximately 180 mm and a MAR165 (at BioCAT) or a 
sample to detector distance of 300 mm and a MAR300 detector (GMCA CAT). XFM 
images were collected at beam line 2-ID-E at the APS as previously described 27 using 
10 keV incident energy, 300 nm step size and 50 ms dwell time per pixel in fly scan 
mode (continuous scan motion in sample x direction). Fluorescence spectra per pixel 
were detected by a four element Si drift detector (Vortex-ME4 by SII 
NanoTechnology) . 

BCDI data sets were collected at beam line 34-ID-C at the Advanced Photon Source 
using instrumentation as previously described 8 . Briefly, data collection involved 
scanning the sample past the x-ray beam while monitoring for the presence of Bragg 
peaks with a relatively low (—0.75 m) sample to detector distance. Once observed, the 
sample to detector distance was usually increased to either 1.0 or 1.8 m for data 
collection at high point-to-point resolution. Crystal size was assessed by scanning 
across the sample in two directions. Full data sets were collected only from crystals 
that were completely irradiated by the incident beam. The data collection strategy 
focused on crystals oriented so that their long axis was approximately parallel to the 
incident x-ray beam making possible reconstruction of crystals with long axes sub- 
stantially longer than the diameter of the incident beam. 

The detector is able to sample a reciprocal space interval of Aq^ = 2n Ap^/XL, 
where Ap^ is the pixel size of the detector, L is the sample to detector distance and X is 
the wavelength of the incident x-ray beam. For the data sets used for image recon- 
struction the pixel size of the detector was 22.5 um, the crystal to detector distances 
were 0.75 m, 1.0 m and 1.8 m and the x-ray wavelength varied from 1.392 A to 1.409 
A (see Supplementary Information, S2). With these parameters the reciprocal space 
samplings, Aq xy , for data recorded on the detector were in the range 0.0056- 
0.0010 nm" 1 . All data sets were collected from a reflection characterized by a scat- 
tering angle of —23° (corresponding to a spacing of —3.65 A) and this reflection was 
scanned by rotating the sample in increments ranging from 0.0017° to 0.0047°. By 
selecting an appropriate value for the rotation increment, the reciprocal space sam- 
pling, Aq z , between consecutive, almost parallel images may be approximately 
matched to the reciprocal space sampling in the detector plane, Aq^. 

The reconstruction process was initiated by associating experimentally measured 
amplitudes with randomly selected phase angles and computing the Fourier trans- 
form. The real- space modification step employed both error reduction (ER) and 
hybrid input-output (HIO) algorithms 1415 for changing the density outside the sup- 
port volume. The initial support volume was generated by smoothing the autocor- 
relation function of the object (i.e. the Fourier transform of the intensities) and taking 
grid points where this function exceeded 4% of the function maximum as defining the 
support 20 . The full image reconstruction protocol required several macrocycles with 
each macrocycle including 10 iterations of ER, 50 iterations of HIO and 50-100 
further iterations of ER. After each macrocycle the support volume was updated using 
the shrink wrap procedure 20 , which allows the support volume to evolve towards the 
shape and size of object that is being reconstructed. Convergence of this process was 
measured by a % 2 statistic that compared the amplitudes computed from the Fourier 
back- transform of real-space arrays modified by the ER or HIO methods with 
experimentally determined amplitudes 6 . After 5-10 macrocycles stable values for the 
support volumes and % 2 statistics were obtained (see Supplementary Information, S3). 

The uniqueness of these solutions to the image reconstruction problem was 
checked by performing replicate runs for each data set in which different sets of 
random phases were used to initiate each run. Compensation was made for the known 
ambiguities, involving translations and inversions of the reconstruction solutions 15 , 
when comparing parallel reconstruction runs. These tests initially showed that the 
shapes and sizes of the crystal reconstructions were consistently determined from the 
absolute values of the final electron densities but the maps of the phase angle com- 
ponent of the complex electron density functions were not necessarily unique at this 
level of data accuracy. In order to obtain more accurate input data for the recon- 
struction process, the values of neighboring intensity measurements were averaged 
together, i.e. the reciprocal- space data sampling was slightly reduced in exchange for 
the greater statistical accuracy of the merged measurements. Reconstructions per- 
formed with this preprocessed data were characterized by lower final % 2 values and a 
much more consistent appearance in the maps of the phase angles associated with the 
complex electron density. Calculations from simulated scattering patterns from 
various objects in which the scattering center was miss-set or in which random 
counting noise was introduced into the intensity distribution were performed to 
assess the impact of these types of error on the resulting complex electron density 
maps. In neither case did the resulting phase maps contain features that resembled the 
distinctive patterns of phase domain separation that are observed in the recon- 
structions from these experimental BCDI data. Also, in test calculations where the 
starting electron density distributions (calculated from random phases) were set to 
purely real values, the expected complex electron density distributions still emerged 
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from the reconstruction process, further demonstrating that the complex electron 
density distributions we have obtained are an inherent consequence of the asym- 
metric intensity distributions observed around the Bragg scattering positions. 

The BCDI reconstructions were initially calculated on a 3D grid defined by two 
dimensions in the plane of the detector (typically 200-300 pixels in side length) and a 
third dimension defined by the sample rocking angle. The resulting maps of the 
complex electron density were subsequently transformed to a cubic Cartesian 
coordinate system for display 8 . 

1. Himmel, M. E. Biomass Recalcitrance: Deconstructing the Plant cell Wall for 
Bioenergy (Blackwell Publishing, Oxford, 2008). 

2. Bradbury, J. H. & Nixon, R. W. The acridity of raphides from the edible aroids. 
Journal of the Science of Food and Agriculture 76, 608-616 (1998). 

3. Prychid, C. J. & Rudall, P. J. Calcium Oxalate Crystals in Monocotyledons: A 
Review of their Structure and Systematics. Annals of Botany 84, 725-739 (1999). 

4. Franceschi, V. R. & Nakata, P. A. Calcium oxalate in plants: Formation and 
Function. Ann. Rev. Plant Biol. 56, 41-71 (2005). 

5. Robinson, I. K., Vartanyants, I. A., Williams, G. J., Pfeiffer, M. A. & Pitney, J. A. 
Reconstruction of the Shapes of Gold Nanocrystals Using Coherent X-ray 
Diffraction. Phys. Rev. Lett. 87, 195505-1-195505-4 (2001). 

6. Williams, G. J., Pfeifer, M. A., Vartanyants, I. A. & Robinson, I. K. Three- 
dimensional Imaging of Microstructure in Au Nanocrystals. Phys. Rev. Lett. 90, 
175501-1-175501-4 (2003). 

7. Robinson, I. & Harder, R. Coherent X-ray diffraction imaging of strain at the 
nanoscale. Nature Materials. 8, 291-298 (2009). 

8. Pfeiffer, M. A., Williams, G. J., Vartanyants, I. A., Harder, R. & Robinson, I. K. 
Three-dimensional mapping of a deformation field inside a nanocrystal. Nature 
442, 63-66 (2006). 

9. Fernandes, A. N. et al. Nanostructure of cellulose microfibrils in spruce wood. 
Proc. Natl. Acad. Sci. (USA) 108, E1195-E1203 (2011). 

10. Tazzoli, V. & Domeneghetti, C. The crystal structures of whewellite and 
weddellite: re-examination and comparison. American Minerologist. 65, 327-334 
(1980). 

1 1 . Frey- Wessling, A. Crystallography of the Two Hydrates of Crystalline Calcium 
Oxalate in Plants. American Journal of Botany 68, 130-141 (1981). 

12. Huang, X., Harder, R., Leake, S., Clark, J. & Robinson, I. Three-dimensional Bragg 
coherent diffraction imaging of an extended ZnO crystal. /. Appl. Cryst. 45, 
778-784 (2012). 

13. Gerchberg, R. W. & Saxton, W. O. A practical algorithm for the determination of 
the phase from image and diffraction plane pictures. Optik 35, 237-246 (1972). 

14. Fienup, J. R. Phase retrieval algorithms: a comparison. Applied Optics. 21, 
2758-2769 (1982). 

15. Fienup, J. R. Reconstruction of a complex-valued object from the modulus of its 
Fourier Transform using a support constraint. /. Opt. Soc. Am. A. 4, 118-123 
(1987). 

16. Millane, R. P. Multidimensional phase problems. /. Opt. Soc. Am. A13, 725-734 
(1996). 

17. Millane, R. P. & Stroud, W. J. Reconstructing symmetric images from their 
undersampled Fourier intensities. /. Opt. Soc. Am. A14, 568-579 (1997). 

18. Miao, J., Charalambous, P., Kirz, J. & Sayre, D. Extending the methodology of Xray 
crystallography to allow imaging of micrometre- sized non- crystalline specimens. 
Nature 400, 342-344 (1999). 

19. Miao, J. & Sayre, D. On possible extensions of X-ray crystallography through 
diffraction-pattern oversampling. Acta Cryst. A 56, 596-605 (2000). 

20. Marchesini, S. et al. X-ray image reconstruction from a diffraction pattern alone. 
Phys. Rev. B. 68, 14010-1-14010-5 (2003). 

21. Rodriguez, J. A., Xu, R., Chen, C.-C, Zou, Y. & Miao, J. Oversampling smoothness: 
an effective algorithm for phase retrieval of noisy diffraction intensities. 

/. Appl Cryst. 46, 1-7 (2013). 

22. Aranda, M. A. G. et al. Coherent X-ray diffraction investigation of twinned 
microcrystals. /. Synchrotron Rad. 17, 751-760 (2010). 



23. Minkevich, A. A., Gailhanou, M., Micha, J.-S., Charlet, B., Chamard, V. & 
Thomas, O. Inversion of the Diffraction Pattern from an Inhomogenously 
Strained Crystal using an Interative Algorithm. Phys. Rev B 76, 104106 (2007). 

24. Petrova, E. V., Gvozdev, N. V. & Rashkovich, L. N. Growth and dissolution of 
calcium oxalate monohydrate (COM) crystals. 

/. of Optoelectronics and Advanced Materials. 6, 261-268 (2004). 

25. Webb, M. A. Cell-Mediated Crystallization of Calcium Oxalate in Plants. 
The Plant Cell. 11,751-761 (1999). 

26. Matsushima, U. et al. Calcium oxalate crystal distribution in rose peduncles: 
Noninvasive analysis by synchrotron X-ray micro-tomography. 
Postharvest Biology and Technology 72, 27-34 (2012). 

27. Twining, B. S. et al. Quantifying Trace Elements in Individual Aquatic Protist 
Cells with a Synchrotron X-ray Fluorescence Microprobe. Anal. Chem. 75, 
3806-3816 (2003). 

28. Whewellite Mineral Data. http://webmineral.com/data/Whewellite.shtml (15 
April 2013) 

29. Pettersen, E. F. et al. UCSF Chimera -- a visualization system for exploratory 
research and analysis. /. Comput. Chem. 13, 1605-1612 (2004). 

30. Schneider, C. A., Rasband, W. S. & Eliceiri, K. W. NIH Image to ImageJ: 25 years of 
image analysis. Nature Methods. 9, 671-675 (2012). 



Acknowledgements 

This work was supported as part of the Center for Direct Catalytic Conversion of Biomass to 
Biofuels, an Energy Frontier Research Center funded by the U.S. Department of Energy, 
Office of Science, Basic Energy Sciences under Award # DE-SC0000997. 1.R. was supported 
by ERC Advanced Grant 227711 'nanosculpture'. C3Bio is supported by the US 
Department of Energy, Office of Basic Energy Sciences. We would like to thank Robert F. 
Fischetti and David Gore for assistance in collection of the WAXS data. GM/CA CAT has 
been funded in whole or in part with Federal funds from the National Cancer Institute 
(Yl-CO-1020) and the National Institute of General Medical Science (Yl-GM-1104). 
BioCAT is a National Institutes of Health-supported Research Center (RR-08630). Use of 
the Advanced Photon Source, an Office of Science User Facility operated for the U.S. 
Department of Energy (DOE) Office of Science by Argonne National Laboratory, was 
supported by the US. DOE under Contract No. DE-AC02-06CH11357. 

Author contributions 

J.B. carried out the final computational analyses and interpretations and wrote the 
manuscript; J.L. collected most of the data and carried out extensive preliminary data 
analyses; R.H. designed the experiments, was involved in all the data collection and advised 
of all aspects of the project; H.I. was involved in analyses of the WAXS data and 
interpretation of the BCDI data; S.C.G. and S.V. collected and analyzed the XFM data; I.R. 
designed the experimental apparatus, was involved in data collection and advised on other 
aspects of the project; L.M. conceived the project, coordinated all aspects and wrote the 
manuscript. 

Additional information 

Supplementary information accompanies this paper at http://www.nature.com/ 
scientificreports 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Badger, J. et al. Three-dimensional Imaging of Crystalline 
Inclusions Embedded in Intact Maize Stalks. Sci. Rep. 3, 2843; DOI:10.1038/srep02843 
(2013). 

0 I This work is licensed under a Creative Commons Attribution 3.0 Unported license. 
I£^^KZ^H To view a copy of this license, visit http://creativecommons.Org/licenses/by/3.0 



SCIENTIFICREPORTS | 3 : 2843 | DOI: 1 0.1 038/srep02843 



5 



